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Abstract 

In the past, a maximum-entropy model was introduced and applied to the study of statistical 
scattering by chaotic cavities, when short paths may play an important role in the scattering 
process. In particular, the validity of the model was investigated in relation with the statistical 
properties of the conductance in open chaotic cavities. In this article we investigate further the 
validity of the maximum-entropy model, by comparing the theoretical predictions with the results 
of computer simulations, in which the Schrodinger equation is solved numerically inside the cavity 
for one and two open channels in the leads; we analyze, in addition to the conductance, the 
zero-frequency limit of the shot-noise power spectrum. We also obtain theoretical results for 
the ensemble average of this last quantity, for the orthogonal and unitary cases of the circular 
ensemble and an arbitrary number of channels. Generally speaking, the agreement between theory 
and numerics is good. In some of the cavities that we study, short paths consist of whispering 
gallery modes, which were excluded in previous studies. These cavities turn out to be all the 
more interesting, as it is in relation with them that we found certain systematic discrepancies in 
the comparison with theory. We give evidence that it is the lack of stationarity inside the energy 
interval that is analyzed, and hence the lack of ergodicity -a property assumed in the maximum- 
entropy model- that gives rise to the discrepancies. Indeed, the agreement between theory and 
numerical simulations is improved when the energy interval is reduced to a point and the statistics is 
then collected over an ensemble obtained by varying the position of an obstacle inside the cavity. It 
thus appears that the maximum-entropy model is valid beyond the domain where it was originally 
derived. An understanding of this situation is still lacking at the present moment. 

PACS numbers: 73.23.-b, 73.63.Kv, 72.70.+m 
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I. INTRODUCTION 



The statistical scattering of waves through open chaotic cavities has been of great interest 
to many groups along the yearsi*^. The investigations that have been carried out are relevant 
to a variety of problems, like the electronic transport through ballistic quantum dots, or the 
scattering of classical waves (e.g., electromagnetic or elastic waves) in chaotic billiards. 

The approach provided by Random-Matrix Theory has been particularly fruitful in the 
study of the statistical fluctuations of transmission and reflection of waves by a number of 
systems, including billiards with a chaotic classical dynamics. Within this approach we wish 
to focus our attention on the model of Refs. POQ, which was introduced originally in the 
context of Nuclear Physics and was then applied to the domain of chaotic cavities. 

We recall that, very generally, we can describe a scattering process in terms of a scat- 
tering matrix S. In the model referred to above, the statistical features of the problem are 
represented by a measure in S-matrix space which, through the assumption of "ergodicity" , 
gives the probability of finding S* in a given volume element as the energy E changes and 
S wanders through that space. The problem is, of course, to find that measure. The key 
assumption is made that in the scattering process two distinct time scales occur, associated, 
respectively, with a prompt, or direct, response due to the presence of short paths, and a 
delayed, or equilibrated, response due to very long paths. It turns out that the prompt, or 
direct, processes can be expressed in terms of the energy average of S, S, also known as the 
optical S matrix. The statistical distribution of the scattering matrix S is then constructed 
through a maximum-entropy "ansatz", assuming that it depends parametrically solely on 
the optical matrix. The notion of ergodicity, which allows replacing energy averages by 
ensemble averages, e.g., {S) = S, is essential to the argument. 

The statistical properties of the conductance predicted by the maximum-entropy model 
we just described have been studied in the past; these predictions have been also compared 
with the results of computer simulations which consist in solving the scalar Schrodinger 
equation numerically for a number of structures^iM. Although in those structures the two 
time scales referred to above were not as well separated as in Nuclear Physics problems, they 
seemed to us to be sufficiently distinct to allow a meaningful description. It is the purpose 
of the present article to investigate further the validity of the maximum-entropy model, by 
extending our earlier studies in the following three ways. 



3 



First, we wish to provide further predictions of our approach for other physical quantities 
in addition to the conductance. For this purpose we analyze the zero-frequency limit of 
the shot-noise power spectrum P at zero temperature. For one open channel (A^ = 1) we 
show that the problem can be reduced to quadratures and, in a number of cases, we can 
even study analytically the influence of direct processes on the average, (P), of the zero- 
frequency shot-noise power spectrum over an ensemble of cavities. For an arbitrary number 
of channels, on the other hand, we show that (P) can be evaluated analytically when direct 
processes are absent [{S) = 0]. 

Second, we wish to extend the computer simulations mentioned above in a number of 
ways: 

i) In some of the cavities used in the present paper the short paths consist of whispering 
gallery modes (WGM), which were excluded in Refs. jsl^ by the type of cavities that were 
used and the way the leads were attached. It is their effecli^ that we wish to describe in 
terms of the optical S matrix which, as we said, is precisely a measure of the short-time 
processes occurring in the scattering problem. Information on the time scales involved could 
be provided by an analysis of the structure of S{E) in the complex-energy plane. Although 
we do not have direct access to the poles of the S matrix, the complex eigenvalues of the 
so-called "effective Hamiltonian" (which essentially consists of the Hamiltonian of the closed 
cavity plus the coupling to the continuum) give evidence of a "sea" of fine-structure, long- 
lived, resonances, plus a collection of shorter lived, more widely separated states. This 
evidence is indicated in the present paper and studied in detail in Refs. 

ii) Earlier numerical simulations were performed for cavities with an applied magnetic field 
(the unitary universality class characterized by the Dyson parameter /3 = 2), in the presence 
of direct processes and for one channel (A^ = 1). The present simulations are performed 
for cavities with time-reversal invariance (the orthogonal universality class, characterized by 
the Dyson parameter /5 = 1), also in the presence of direct processes and for one (A^ = 1) 
and two (A^ = 2) open channels. 

Third, we shall pay closer attention to the discrepancies between theory and numerical 
experiments. Indeed, discrepancies similar to the ones that we shall observe in this paper 
were already present, to a certain extent, in Ref. P|, but were overlooked at that time. 

The paper is organized as follows. In the next section we first give a brief presentation 
of the maximum-entropy model, recalling the assumptions that are used in its derivation; 



4 



these considerations will be important in the discussion to be presented in Sec. IIVI We 
then study a number of predictions of the model with regards to the statistical properties 
of the conductance and the shot-noise power spectrum at zero temperature. In Sec. 11111 
we present the results of the numerical simulations and the comparison with theory. Sec. 
1111 Al is devoted to the one-channel case (A^ = 1) and Sec. 1111 Bl to two channels (A^ = 2). 
Finally, we discuss our results in Sec. IIVI putting particular emphasis on the discrepancies 
found between theory and numerical simulations. We include an appendix, where some of 
the algebraic details of the relevant one- and two-channel statistical distributions are given. 



II. STATISTICAL MODEL FOR THE DESCRIPTION OF QUANTUM CHAOTIC 
SCATTERING IN BILLIARDS 

We present below the main ideas behind the maximum-entropy model briefly described 
in the Introduction. This model was introduced in the past in the domain of Nuclear 
Physics and was later used to study the quantum mechanical scattering occurring inside 
ballistic cavities (whose classical dynamics is chaotic) connected to the outside by means of 
waveguides^>^. 

The scattering problem can be described in terms of a scattering matrix S. If the cavity 
is connected to two waveguides supporting A^ channels each, the dimensionality of the S 
matrix is 2N. As we mentioned in Sec. HJ the model proposes a measure in S'-matrix space 
which, through the assumption of ergodicity, describes the probability of finding in a given 
volume element as the energy E changes and S wanders through that space. We write such 
a probability as 

dPl^,^^{S)=pf,\{S)df,p{S), (2.1) 

where p|^j(5'), referred to as the probability density, depends parametrically on the optical 
matrix {S), as detailed below. In the above equation, dn^^\S) is the invariant measure 
for the universality class (3 [we shall assume throughout that J d^p{S) = 1]. Here we shall 
consider the cases (3 = 1 (the orthogonal case) and (3 = 2 (the unitary case), corresponding 
to cavities with and without time-reversal invariance, respectively, and in the absence of 



spin. The problem is, of course, to find 



mm 



1^) 



[S). To this end, a number of assumptions are 



made, as we now explain (see Refs. UOUIH). 

1) The study of the statistical properties of S{E) over an ensemble of cavities is simplified 



by idealizing S{E), for real E, as a stationary random (matrix) function of E satisfying the 
condition of ergodicity. 

2) As explained in Sec. El we assume that our scattering problem can be characterized in 
terms of two time scales, arising from the prompt and equilibrated components; the prompt 
response can be described in terms of the averaged S matrix (S), also known as the optical 
S matrix. 

3) We assume E to be far from thresholds, so that, locally, S{E) is a meromorphic 
matrix function which is analytic in the upper half of the complex-energy plane and has 
resonance poles in the lower half plane. From this follow what we have called in the past 
the "analitycity-ergodicity" (AE) properties: 



This expression involves, on its left-hand side, only S matrix elements, whereas S* matrix 
elements are absent; on the right-hand side, only the optical matrix (S) appears. More 
generally, if f{S) is a function that can be expanded as a series of non- negative powers of 
the S matrix elements, we must have the reproducing property^ 



is such that the average S matrix is the optical matrix (S), the AE requirements ()2.2|) and 
hence the reproducing property ()2.3p are satisfied, and the entropy S[p] associated with it, 
S[p] = — f P(s){S) In p(^s){S)dfi{S), is greater than or equal to that of any other probability 
density satisfying the AE requirements for the same (S). 

With regards to the information-theoretic content of Poisson's kernel, we have to distin- 
guish between i) general properties, like unitarity of the S matrix (flux conservation), ana- 
lyticity of S{E) implied by causality, and the presence or absence of time- reversal invariance 
(and spin-rotation symmetry when spin is taken into account) which determines the univer- 
sality class (orthogonal, unitary or symplectic), and ii) particular properties, parametrized 
by the ensemble average (S), which controls the presence of short-time processes. System- 
specific details other than the optical matrix (S) are assumed to be irrelevant. The optical 
matrix {S) is the only "physically relevant parameter" assumed in the model. 




(2.2) 



(/(^)) = /((^))- 



(2.3) 




(2.4) 
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From the probability distribution of Eqs. ()2.ip and ()2.4p one can find the statistical 
properties of the quantities of interest over an ensemble of cavities. In this paper we shall 
be concerned with the conductance and the zero-frequency shot noise power spectrum. 

The dimensionless dc conductance [g = G/ {e^/h)] at zero temperature and for the spinless 
case is given by Landauer's formula!^ 

g = T = tr{tt^) 

= E^- (2.5) 

a 

where Ta (a = 1, ■ ■ ■ , N) are the eigenvalues of the Hermitean matrix tt^ , and the transmis- 
sion matrix t is an x block of the 2A^-dimensional S matrix which, in turn, is written 

as 

r t' 

S= . (2.6) 

t r' 

The zero-frequency limit of the shot-noise power spectrum at zero temperature can be 
expressed aaM^ 

2p2 

P = PoJ2^a{l-ra), Po = 2eV—. (2.7) 

a=l 

The average of P over an ensemble of cavities will be written in the two alternative ways: 
(P) = Po(Er,(l-r,)\ (2.8a) 



, n=l 



{Pp) ' (Pp) = ^eV—{T). (2.8b) 



Here, Pp is the result that would obtain if the noise were a Poissonian process, i.e., if there 
were no correlations among electrons and the electronic transport were completely random; 
T is the dimensionless conductance, Eq. ()2.5|1 . We see that since the shot-noise power is not 
determined simply by the conductance, it is only in the limit ^ 1 (a = 1, ■ ■ ■ , A^) that 
we recover the Poissonian result. 

It is clear that we need, for our purposes, the joint probability distribution of the r^'s. 
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This can be found from Eq. ()2.4|) as 

w;jj}(ri, . . . ,rAr) = Ci 



(1) / ^ \ _ ^a<b - n\ 




w}^\ {ti,...,tn) = C2Y\_ i^a - n) 

a<b 



det(/- (g)(g)t)|^+i 
|det(/- (5)5t)|2JV+i 

det(/- (g)(g)t)p^ 
|det(/- (5)5t)|4^ 



i=l 



(2.9a) 



(2.9b) 



for (3 = 1 and (3 = 2, respectively. The quantity is a normahzation constant. The 
unitary matrices w^*-* are the ones that occur in the polar decomposition of the S matrij*^ 



S 









yi3) 

v^'^ 



(2.10) 



where r stands for the N x N diagonal matrix constructed from the the eigenvalues 
(a = 1, ■ ■ ■ , A^) of the Hermitian matrix tf^ [see Eq. fl2.6p ] and the w*^*^ are arbitrary N x N 
unitary matrices for P = 2, with the restrictions v^^^ = [v^-^^]^ and v^^^ = [v^'^^]'^ for P = 1. 

In what follows we study, in particular, the cases in which the two waveguides connecting 
the cavity to the outside may support one, two, or an arbitrary number of open channels. 



A. The one-channel case, A'^ = 1 



In this case we have only one r, which coincides with the conductance T, whose probability 
distribution can thus be written from Eqs. ()2.9|) as 



, (1) _ 1 /'^ /'^ |det(/-(g) (g)t)|3/2 dadp 
^^^^^^-Vfio L |det(/-(5)5t)|3 (2^' ^2.11a) 



_ r r Idet(/-(g) {Syr dad(3d^dS 
^'^^ ^ - X "Jo |det(/-(5)5t)|4 (2vr)^ " ^^"^^^^ 

The polar representation of S for = 1 is written down explicitly in Eq. (jAlll of the 
Appendix. 

In the absence of direct processes, i.e., {S) = 0, the T distribution of Eqs. ()2.1ip reduces 
to the well known results 

A'\t) = (2.12a) 
w^^\t) = 1, (2.12b) 
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for the orthogonal (/3 = 1) and unitary (/3 = 2) cases, respectively. 

The T distribution for the unitary case, Eq. ()2.11b|) . can be integrated explicitly^. As 
an example, for the particular case (t) = {t') = 0, corresponding to direct reflection and 
no direct transmission, and assuming, for simplicity, the "equivalent-channel" case {\{r)\ = 
l(r')l), one finds 

[(1- |(r)|2)2 + 4|(r)|2T] / 

For the case of direct transmission and no direct reflection, the result is obtained from the 
previous equation by replacing | (r) | by | (t) \ and T by 1 — T. 

The r(= T) distribution for the unitary case given in Eq. ()2.1H|1 allows us to study the 
effect of direct processes on the averaged shot-noise power spectrum (P) of Eq. ()2.8b|) : 
this case is particularly suited to gain some physical insight, since the result for (P) can 
be expressed analytically in a remarkably simple fashion. For the particular case of direct 
reflection and no direct transmission ((t) = {t') = 0), and assuming |(r)| = |(r')|, one finds, 
from Eq. ()2.1Hj) . the result: 

(P)(^) _ 1 5-9|(r)|^ + 4|(r)r . . 

(Pp)(2) 5 3-4|(r)|2 + |(r)|4 ' ^ ' ' 

while for direct transmission and no direct reflection ((r) = (r') = 0), and assuming \ {t) \ = 
I {t') I , one obtains 

(P)(2) 1 5 - 9\{t)\^ + mf ^2.15) 



(Pp)(2) 5 3 + 4|(t)P-|(t)|4- 
In Fig. ^ the behavior of the ratio {P)^'^^ {Pp)^'^'^ as a function of |(r)| = \{r')\ for the 
direct reflection case {{t) = {t') = 0), Eq. ()2.14|) . is shown as the upper solid curve; the 
lower solid curve shows the case of direct transmission as a function of = \{t')\ (when 
(r) = (r') = 0), Eq. (j2.15|) . For the upper curve, the ratio {P)^'^^ {Pp)^"^^ increases as a 
function of \{r)\] since, as |(r)| 1, w{T) S{T), at first sight one would expect, in this 
limit, the ratio {P)/{Pp)^'^^ to increase towards the Poissonian value unity. That this is not 
the case is due to the fact that both (T) and (T^) tend to zero linearly with 1 — |(r)| as this 
quantity tends to zero. 

For the orthogonal symmetry class (/3 = 1) we have not succeeded in finding an analytical 
expression for the conductance distribution, even for the particular cases studied above. For 
these cases, the ratio (P)*^^Y(Pp)*^^^ was thus calculated numerically from Eq. ()2.11a|l and 
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the results are also presented in Fig. for comparison with the unitary case; we observe 
that the ratio (P)*-'^-*/ (Pp)'-^-' is always larger for (3 = 1 than for (3 = 2. 

We wish to point out a property of the average shot-noise power (P)*^^) of Eq. (j2.8ap . in 
the present one-channel case. Poisson's kernel of Eq. (j2.4j) has the property that has been 
called "covariance"— : p|^|(S') = where S = UqSVq, Uq and Vq being fixed unitary 

matrices for (3 = 2, with Vq = t/^f for /3 = 1, the same transformation being applied to the 
optical (S). The invariant measure is invariant under this transformation. For /5 = 2, one 

1 



can verify that the unitary matrices Uq 



1 



and Vo = 1 exchange r (r') and t {t') and 



their corresponding average values appearing in (S). For the case (3 = 1, we have t = t'. 
If we also have r = r', as in the case of a system with "left-right symmetry", the matrix 

1+i 1 -i 



1 

V2 



1-i 1+i 
above transformations 



switches r and t and the corresponding optical parameters. The 



ieeps P = Pot(1 — r) invariant. As a consequence, (P)^^^ remains 
invariant under the interchange (r) ^ (t), (r') ^ {t') for (3 = 2, and (r) ^ (t) for the 
particular (3 = 1 case mentioned above. We observe that, indeed, the numerators of Eqs. 
fl2.14|l and ()2.15|1 . which are proportional to {P)^'^\ do fulfill this property. However, for the 
(3=1 case considered here, this symmetry does not apply. 

In the present one-channel case one can write down an expression for the distribution of 
the "dimensionless" shot-noise power spectrum [see Eq. ()2.7|) ] 



V 



r 



(2.16) 



which lies in the range < rj < ^ (we are using the notation of Ref. Il3f ). Since is a 
function of the conductance, we can make an elementary change of variables and write 

■^(/3)(r) 



1 

T = - 

2 



|1 2r|j^^^(^) 
1 ± v^l -4r/ 



(2.17a) 
(2.17b) 



Thus the distribution in question is given by: 



+ w 



2 ) 



(2.18) 



where w'^^^{t) is given in Eqs. ()2.1H) . For (S*) = 0, the result of this last equation ()2.18|) 
reduces to Eq. (95) of Ref. |l3. 
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FIG. 1: The ratio {P)/{Pp) as a function direct reflection |(r)| = \{r')\ (indicated in the upper 
horizontal line as the abscissa) for the case {t) = {t') = is shown as the two upper curves. The 
two lower curves show the same ratio as a function direct transmission \ {t)\ = \ {t')\ (indicated in 
the lower horizontal line as the abscissa) for the case (r) = {r') = 0. The dashed lines correspond 
to the orthogonal universality class {(3 = 1) and the solid lines to the unitary class (/3 = 2). 

B. The two-channel case, N = 2 

In the two-channel case the matrix tt"^ is two-dimensional and has two eigenvalues Ti,T2, 
whose joint probability distribution can be written from Eqs. (12 .yj) as 



^''^^ ^ 4 J J \det{I - S {S)^)\^ J t-K I 

(2.19a) 
(2.19b) 

Here, djji{v^'^^) is the invariant measure for the unitary matrices used to represent 5* in 
its polar form, Eq. ()A3|) : the explicit form of dfi{v^^^) is given in Eqs. ()A6|) and ()A7|) . 

From the above expressions we can evaluate the probability distribution of the conduc- 



11 



tance as 



^\s)iT)= / w\',lin,T - n)dn, (2.20) 



and the ratio (P) / (Pp) for the shot-noise power spectrum as 

iP"}' (EL-.) ■ ' ' 

In the absence of direct processes, {S) = 0, we obtain for Wq^\ti, T2) the well known results^: 

4^)(n,r.) = I^L^ (2.22a) 

wi^\ruT,) = 6{n-T2)\ (2.22b) 

and for the conductance distribution w\^\t) 

m f 1^, < T < 1 

wi^\T) = { ' (2.23a) 

1^ I (T - 2y/T^) , 1 < T < 2. 

wi^\T) = 2[1 - \1-T\f. (2.23b) 



C. The case of arbitrary N 

In the absence of direct processes, {S) = 0, various results concerning the average and 
variance of the conductance are known^ and will not be reproduced here. 

Not known, to our knowledge, is the behavior of the shot-noise power spectrum for 
arbitrary A^, even for {S) = 0. We calculate below, for such a situation, the average (P) for 
the orthogonal and the unitary cases. 

The numerator of ()2.8b|) can be written as 

N 

- T^ \ = {tr{tt^)) - {tr{ttht^))^ 

N N 

X^d^afeHo^^- ^ {tabtcdt*cbtla)(f^ 
a, 6=1 a,b,c,d=l 

N N 

a, 6=1 a,b,c,d=l 

12 



, a=l 




The notation (■ ■ ■ )o^'' indicates an average over the invariant measure for the universahty 
class (5. 

In the last line of Eq. ()2.24|) the upper indices 21 indicate the 21 block of the S matrix 
in Eq. (fTK|l . 

Averages of monomials of the type 



(2.25) 



10 and [3, for /3 = 1 and (3 = 2, respectively. We now consider these 



were studied in Ref. 
two cases separately. 

In the orthogonal case, /3 = 1, we denote Q{1) = M, just as in Ref. In that reference 
one finds the results 



0. 



m: 



2iV+ 1 



A 



a/3 



a/3,75 



A 
+B 



M:/Mi' + M^jM;'/ 



a'5' ..P'l' 



0/3 75 



where 



A 



{2N + 1){2N + 2) 



B 



Q/3 75 



2A^+ 1 



1/3 75 



(2.26a) 



(2.26b) 



(2.26c) 



2iV(2A^ + 3) ' 2iV(2A^ + 3)' 

Substituting the results ()2.26|1 in Eq. (|2.24j) we find for the average of P, Eq. ()2.8b|l . for 
the orthogonal case: 

(iV+l)2 



{P) 



(/3=1) 



{PP) 



(/3=1) 



2iV(2iV + 3) 

In the unitary case, /? = 2, we denote (5(2) = Q, just as in Ref. 



one finds the results 



Q 



Ql 



al3 
■/3,75 



iV 



(2.27) 
In that reference 

(2.28a) 



1 



{2Ny - 1 L 
1 



Al^'A^'/' + AI'^'A^'^' 



^07 "/35 



^07 '^/3(5 



07 /3(5 ' a.'y /3o 



2N[{2Ny - 1] 

which has to be substituted in Eq. ^TI^ . For (P)(^=2), Eq. we find: 



(2.28b) 



4Ar2 _ I 



(Pp)(^=2). 



(2.29) 
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For a large number of open channels, ^ 1, Eqs. (j'2.'27p and (j2.29|) give {P)^^^ 



3Q 



KPp)^^) ^ NPo/8, just as in Refs. 

The ratio {P)^'^'^/ {Pp)^f^^ from Eqs. (IT?7|) and (IT^ is plotted in Fig. Has a function 
of the number of channels A^. We observe that this ratio is always larger for the orthogonal 
{(5 = 1) than for the unitary case (/? = 2), just as was noticed in the results shown in Fig. 
^for the one-channel case. This effect indicates that time reversal symmetry pushes the Tq 
distribution towards small r^'s [for = 1 this effect is given by Eq. (|2.12ap ] in such a way 
that (P)(^=^) gets closer to Poisson's value. 



0.4 



V 

A 0.3 



0.25- 



1 ' 

1 


1 1 1 
-|3=1 




- P = 2 






1.1,1,1 



2^ 

2 4 6 8 10 

N 



FIG. 2: The ratio {P) / {Pp) for (3=1 (upper curve) and [3 = 2 (lower curve) for {S) = 0, as a 
function of the number of channels A^. 



III. NUMERICAL SIMULATIONS 

The maximum-entropy approach that we have been discussing is expected to be valid 
for cavities in which the classical dynamics is completely chaotic -a property that refers 
to the long-time behavior of the system- as in such structures the long-time response is 
equilibrated and classically ergodic. 

In Refs. the scalar Schrodinger equation was integrated numerically for a number 
of 2D cavities in order to examine to what extent our approach really holds. In those 
references the analysis was performed for the conductance distribution w{T). The cavities 
were subjected to a magnetic field {(3 = 2) and they were connected to the outside by 
waveguides admitting one open channel (A^ = !)• Moreover, the structures were such that 
they obviously supported short paths associated with direct reflection from a barrier, direct 
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transmission from one lead to the other, or skipping-orbit trajectories in the presence of the 
magnetic field. 

In what follows we consider the numerical solution of the Schrodinger equation for 2D 
structures which again support direct processes. Now the system is not immersed in a 
magnetic field, so that it is time- reversal invariant {(3 = 1). We mainly study the one- 
channel case, = 1 (Sec. IIIIAI below), although we also present some results for N = 2 
(Sec. ITTTbI) . 

In addition to the conductance distribution, the average of the zero-frequency shot noise 
power spectrum is also studied, in order to examine further the applicability of the model. 
Ensembles of similar systems are obtained by introducing an obstacle inside the cavity and 
changing its position (see Figs. |21 ID and [7| below). In all cases the optical 5* matrix (S) was 
extracted from the data and used as an input in Eq. ()2.4j) . or in the various results of the 
Sec. |nj to produce the theoretical predictions to be compared with numerical experiments. 
. n this sense all of our fits are "parameter free" . For details of the numerical study see Refs. 

303- 



A. The one-channel case, = 1 

When the energy E lies inside the interval 2mw2 [A^^vr^, (A^ + l)^vr^], each waveguide (of 
width W) supports exactly A^ open channels. In units such that h'^/2mW'^ = 1, one open 
channel (A^ = 1) occurs for E e [7r^,47r^] ^ [10,40]. We need to study S{E) in energy 
intervals AE' not too close to either threshold, in order to avoid threshold singularities. 



1. Statistical properties of the conductance 

Fig. 121 shows, as insets, the structures for which the numerical study was performed: 
they consist of a Bunimovich stadium connected to two waveguides directly, as in panels 
(a), (b) and (c), or through a smaller half stadium, as in (d). The structures are spatially 
asymmetric. 



The histograms were obtained by solving the Schrodinger equation inside these structures 
and collecting the data in the energy interval /S.E = [22, 23] (in the units explained above). 
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FIG. 3: The distribution W (not be confused with the width of waveguide W used in the text) of the 
conductance. normaUzcd to the total number of cases, for the structures shown in the insets. Each bin 
shows the frequency that occurred in that interval. The histograms were obtained from a numerical solution 
of the Schrodinger equation in the energy interval E G [22,23], as explained in the text. The theoretical 
distributions, obtained from Poisson's kernel using the optical S matrix extracted in each of the four cases, 
were calculated at 15 points (in the interval < T < 1) which were then joined to obtain the continuous 
curves. The agreement between theory and numerical simulations is, in general, good; the largest deviations 
occur in panel (d), where the optical, or direct, transmission, is the largest, due to whispering gallery modes 
in the small cavity. 

and then across an ensemble of 200 positions of the obstacle, which is also shown in the figure. 
In that energy interval, 20 equally-spaced points were considered: these points are farther 
apart than the correlation energy, as it appears from the negligible correlation coefficient 
(over the ensemble) that was obtained for the transmission and reflection amplitudes for two 
successive points. The optical S matrix, obtained as an energy plus an ensemble average of 
S, i.e., (S), was extracted from the data and the optical reflection and transmission matrix 
elements are given in Table HJ 

The optical S matrix was substituted in Eq. (j2.11aj) . which is the theoretical prediction 
for the conductance probability distribution w^^l{T), giving, after normalizing to the total 
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TABLE I: The optical reflection and transmission matrix elements for the four cases in Fig. |31 



Case 


if) 


(i) = (to 


{?) 


Fig. 3a 


0.0007-0.0651i 


-0.0725+0.0078i 


0.0040+0.0008i 


Fig. 3b 


-0.0384-0.0213i 


-0.0767-0.0211i 


-0.0388-0.0210i 


Fig. 3c 


0. 1462-0. 0242i 


-0.1236+0.3672i 


0.0375+0.1035i 


Fig. 3d 


0.1106-0.1581i 


-0.0591-0.6055i 


0.2586-0.0054i 



number of cases, the results shown as the continuous hues in Fig. 3. In other words, the 
parameters on which the theoretical results depend, i.e., the optical S matrix elements, were 
not obtained by a variation procedure designed to find a best fit, but rather extracted from 
the numerical experiment. In panels (a) and (b) the optical 5* matrix is very close to zero, 
indicating a negligible presence of direct processes, so that the continuous curve in both cases 
is practically given by Eq. ()2.12a|) : we mainly have long lived states in these structures. 
The elements of the optical S matrix grow larger as we proceed to the remaining panels. 

The agreement between theory and numerical experiments is very good for (a) and (b). 
One observes some deviations in panel (c); the deviations are largest for panel (d), where 
the optical, or direct, transmission, is largest. In (c) the direct path between the two waveg- 
uides is obvious. In (d) the direct path occurs inside the smaller stadium, which supports 
whispering gallery modes, while the larger stadium provides a "sea" of fine-structure states. 
This effect is seen in Fig. |3J where a plot of the square of the scattering wave function for a 
fixed energy E, i.e., |'?/'£;(r)p, exhibits a concentration along the wall of the smaller stadium. 
Indeed, the reason for attaching a small stadium to a large one in Fig. Efd) is precisely 
that the whispering gallery modes which have been seen in calculations for small cavities, 
as in Refs. 00], are gradually destroyed if the size of the cavity is increased, because of the 
long way the whispering gallery mode would have to traverse (for more details the reader is 
referred to Ref. Q)- This seems to be the reason why no effects from WGM are seen in the 
stadium in Fig. Efb). 

We wish to investigate the case of Fig. 01 (d) further. For the optical S matrix extracted 
from the data, the probability of small T's predicted by the theory is larger than that 
found in the numerical simulation, and vice versa for T ~ 1. This effect is not a statistical 
fluctuation, but rather a systematic discrepancy, as it was observed in most cases where the 
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FIG. 4: The square of the scattering wave function, i.e., |?/jE(r)p, for a fixed energy. The geometric structure 
consists of a small stadium coupled to a larger one. The geometry is the same as that shown in Fig. Old). 
We interpret the concentration of the wave function along the wall of the small cavity as a whispering gallery 
mode. The larger stadium provides a "sea" of fine-structure resonances. 

transmission part of the optical S matrix is as large as that occurring in Fig. Efd). To find the 
origin of the discrepancy we have to realize that, in order to apply our model meaningfully, 
an energy interval AE over which we do statistics must be such that the "local" optical 
matrix {S{E)) is "reasonably constant" across it, while, at the same time, such an interval 
should contain many fine structure resonances: in other words, the notion of "stationarity" 
should be approximately valid. Figure |S1 shows the ensemble expectation value {t{E)) of 
the transmission amplitude t{E) as a function of the energy E, for the structure shown in 
Fig. El (d). Although this quantity is certainly not expected to be literally constant, the 
question is whether its variation across the energy interval AE = [22, 23] (used in Fig. El 
(d)) can be considered to be the cause of the discrepancy that we have seen between theory 
and numerical experiment: the following results support our belief that the answer to this 
question is positive. Figure IHl shows the conductance distribution for the same structure of 
Fig. El (d), obtained using a number of energy intervals twice as small. For instance, panels 
(c) and (d) correspond to the two subintervals in which the original one, AE = [22, 23], was 
divided. Panels (a) and (b) show the data for two other similar subintervals. We observe 
that the agreement is now very good. 
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FIG. 5: The real (solid line) and imaginary (dashed line) parts of the ensemble expectation value {t{E)) 
(an element of the "local" optical matrix) as a function of energy, for the structure shown in Fig. |31(d). The 
c^uestion is whether the variation of these quantities inside the energy interval AE ~ [22, 23] is responsible 
for the discrepancy seen in Fig. 13 (d) . 

2. Statistical properties of the shot-noise power spectrum 

The shot-noise power spectrum at zero temperature of Eq. ()2.8b|) was calculated, both 
numerically as well as from our theoretical model, for the same structures shown in Figs. El 
and El For one channel, = 1, Eq. (|2.8bjl simplifies to 



so that in this case we do not have more information than that contained in the conductance 
distribution. However, for completeness, we present the results in Table ITTl 

Notice that the results in the first two rows of the table, i.e., those arising from Figs. 
Of a), (b) (with an optical 5* close to zero), compare well with the prediction of Eq. ()2.27|1 
for = 1, i.e., (P)/(Pp) = 0.4, and with Fig. ^for /3 = 1 and (S) = 0. 

Notice also that from row 1 (or 2) to row 4 of the table the optical transmission increases 
and {P)/{Pp) decreases, in qualitative agreement with the result of Fig. Q for = 1, (r) 
literally equal to zero and increasing (t). 

B. The two-channel case, N=2 

In units such that h?/2mW'^ = 1, two open channels {N = 2) occur in the energy interval 
E e [4-7r^, Qvr^] [40, 90]. In view of the experience gained in the one-channel case described 



(P) _ (T(l-T)) 



(3.1) 
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FIG. 6: The distribution W of the conductance, normaUzed to the total number of cases, for the same 
structure as in Fig. 121(d), but using energy intervals, shown in each panel, twice as small for the construction 
of the histograms. Panels (c) and (d) show the same data of Fig. 121(d), but analyzed inside each of the 
two subintervals. Panels (a) and (b) show the data for two other similar subintervals. The agreement with 
theory is very good. 

above, the energy interval AE was literally reduced to a point, and the statistical properties 
of the conductance and the shot-noise power spectrum were studied across the ensemble for 
a fixed energy E: more specifically, 200 samples were collected at E = 75, varying, just as 
before, the position of the obstacle inside the cavity. 

1. Statistical properties of the conductance 

The numerical simulation was done by solving numerically the Schrodinger equation for 
the structures shown in Figs. Efa) andEfd). The theoretical prediction for the conductance 
distribution w{T) was obtained using Eq. ()2.20|) . which in turn makes use of Eq. (j2.19a|) 
and the equations given in the appendix: the integrations occurring in the equations of the 
appendix were performed numerically using a Monte Carlo method (Metropolis algorithm). 
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TABLE II: The shot- noise power spectrum of Eq. H3.1|) . = 1 



Case 


Numerical 


Theoretical 


17i cr '-ia 




4000 






4000 


rig. OL 


n 9070 + n nni 


39^1 


Fie 3d 


2315 + 0026 


0.2438 


Fig. 6a 


0.1765 ±0.0031 


0.1959 


Fig. 6b 


0.1972 ±0.0041 


0.2001 


Fig. 6c 


0.2576 ± 0.0042 


0.2587 


Fig. 6d 


0.2104 ± 0.0029 


0.2187 



The optical S matrix {S{E)) that was used was extracted from the data at i? = 75 and 
across the ensemble; it is given in the following equation for Fig. Ufa) 



{S{E)) = 



-0.0312 + 0.02592 


0.0805 - 


- 0.0393i 


-0.1391 - 


f 0.0599i 


0.1186- 


{- 0.0908i 


0.0805 - 0.0393i 


-0.1600 


± 0.2470i 


-0.0763 


- 0.009U 


0.0982- 


\- 0.0056i 


-0.1391 ±0.0599i 


-0.0763 


- 0.0091i 


0.1032 ± 


-0.0091i 


0.0357- 


0.05562 


0.1186 ±0.0908i 


0.0982 H 


- 0.0056i 


0.0357- 


- 0.0556i 


0.0723 ± 


0.0764i 



(3.2) 

while for the cavity shown in Fig. W[h), the optical S matrix is given by 



{S{E)) = 



0.1536 - 
0.0469 ± 


0.1256i 
0.0313i 


0.0469 ±0.0313i 
0.0777 - 0.0255i 


0.0703 - 
0.0204- 


- 0.4275^ 
f 0.2628i 


0.0743 - 0.2620Z 
0.7589 ±0.1623i 


0.0703 - 
0.0743 - 


0.4275Z 
0.2620i 


0.0204 ±0.2628z 
0.7589 ±0.1623i 


-0.0135- 
0.0929 - 


- 0.1261Z 
0.0125i 


0.0929 -0.0125Z 
-0.0452 - 0.0195i 



(3.3) 

The blocks in the two previous equations indicate the optical transmission and reflection 
matrices, as in Eq. ()A5|) . Notice that {S{E)) ^ in Eq. (j3.2p . while Eq. (j3.3|) shows a large 
direct transmission. 

The results given in Fig. [7| show a strong effect of direct processes on the conductance 
distribution, whose trends are well described by Poisson's kernel. 
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FIG. 7: The distribution W of the conductance, normalized to the total number of cases, for the structures 
shown in the insets and for two open channels [N = 2): the structures in panels (a) and (b) correspond 
to those shown in panels (a) and (d) of Fig. |31 respectively. The structure in (a) consists of a Bunimovich 
stadium connected to two waveguides directly, while in (b) the connection is done through a smaller half 
stadium. The histograms were obtained from a numerical solution of the Schrodinger equation for these 
structures at the energy E = 75, and constructing an ensemble of 200 positions of the obstacle. The optical S 
matrix was extracted from the data and used to obtain, from Poisson's kernel, the theoretical distributions; 
these were computed at 50 points in the interval < T < 2, which were then joined to obtain the continuous 
curves. The trends shown by the numerical distributions are well bescribed by theory. 



2. Statistical properties of the shot-noise power spectrum 

The theoretical predictions for the average of the shot-noise power spectrum of Eq. ()2.21|) 
were compared with the results of the numerical simulation. Notice that in the two-channel 
case the statistics of the shot-noise power spectrum gives information which is not contained 
in the conductance distribution. The comparison is shown in Table lllll for the same cases 
denoted clS (^db ) and (b) in Fig. [7| We note that ), whose (S) is close to zero, compares 

reasonably well with the theoretical result [{P)/{Pp) = 0.31] of Eq. ()2.27p for = 2 and 
{S) = 0. 
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TABLE III: The shot-noisc power spectrum of Eq. (|T?T|) . N = 2 



Case 


Simulation 


Theoretical 


(a) 


0.2959 ± 0.0030 


0.3300 


(b) 


0.1144 ±0.0011 


0.1200 



IV. CONCLUSIONS AND DISCUSSION 

The statistical properties of the dc conductance in chaotic cavities have been investigated 
in the past in the framework of the maximum-entropy model described in the previous 
sections. Within the same framework, in the present paper we have gone further by studying, 
in addition to the conductance, the zero-frequency shot-noise power spectrum. The shot 
noise is a more complicated quantity than the conductance, in the sense that it involves 
electron correlations due to the Pauli principle. We have been particularly interested in the 
effect that direct processes consisting of whispering gallery modes have on the conductance 
and on the shot-noise power; these modes were promoted by choosing properly the structure 
of the cavities and the position of the leads. This kind of direct processes were, in fact, 
avoided in previous publications by some of the present authors. For the two symmetries 
{(3 = 1,2) studied here we have found that the ratio {P)^^^ {Pp)^^\ as a function of the 
number of channels for (5*) = 0, is larger for /5 = 1 than for (3 = 2, indicating that small 
values of the transmission eigenvalues are favored by time-reversal symmetry. 

We have found that the agreement between the theoretical predictions and the results of 
computer simulations performed for one and two open channels is generally good. However, 
the systematic discrepancies that we have observed lead us to revise the notions under which 
our model has been constructed. 

Indeed, the maximum-entropy model described in Sec. |n] relies on a number of assump- 
tions. For instance, the extreme idealization is made of regarding S{E) as a "stationary 
random (matrix) function" of energy. As a consequence, the optical matrix {S{E)) is con- 
stant with energy and the characteristic time associated with direct processes is literally zero. 
The property of stationarity allows defining the notion of ergodicity which, together with 
analitycity, gives the reproducing property, Eq. ()2.3|) . which is essential for the definition of 
Poisson's kernel (PK) of Eq. ^T^i . 

Needless to say, in realistic dynamical problems stationarity is only approximately ful- 
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filled, so that one has to work with energy intervals AE across which the "local" optical 
matrix {S{E)) is approximately constant, while, at the same time, such intervals should con- 
tain many fine-structure resonances. This compromise can actually be realized in Nuclear 
Physics, where the optical (S) arises from the tail of many distant resonances or from a 
single-particle resonance that lies so far away in the complex-E' plane to act as a smooth 
background on top of the fine-structure compound-nucleus resonances: hence the huge sepa- 
ration between the two time scales. In contrast, as we saw in Sec. IIIIl such a compromise is 
difficult to fulfill for the physical systems studied here: this we believe to be the origin of the 
discrepancies observed between theory and numerical simulations. (Indeed, discrepancies 
similar to the ones that we have observed in this paper were already there, to a certain ex- 
tent, in Refs. but were overlooked at that time.) In the present paper we give evidence 
that reducing AE literally to a point and collecting data over an ensemble constructed by 
changing the position of the obstacle inside the cavity, the agreement between theory and 
experiment is significantly improved, being excellent in several cases. In other words, PK 
gives a good description of the statistics of the data taken across the ensemble. 

It is interesting to remark that also in Ref. cases had been found in which stationarity 
obviously did not hold. Energy averages were out of the question in those cases, so that an 
ensemble was generated by adding "noise" along the wall: it was found that PK gave an 
excellent description of the data collected across the ensemble at a fixed E. This point was 
merely indicated at that time and no results were published. 

Thus the results shown in the present paper give evidence that PK is valid beyond 
the situation where it was originally derived, which required the properties of analyticity, 
stationarity and ergodicity, plus a maximum-entropy ansatz. It is as though the reproducing 
property of Eq. ()2.3p were valid even in the absence of stationarity and ergodicity (analyticity 
is always there, of course). Even at the present moment we are unable to give an explanation 
of this fact. A few remarks are in order in connection with this point. 

Let us take the invariant measure dfi{So) of Sec. |n] as a model for the description of 
scattering by a chaotic cavity described by the scattering matrix 5*0 and assumed to have ideal 
coupling to the leads. Brouwer has shown (see Ref. 16[, Sec. V) that when such a chaotic 
cavity is coupled to the leads by means of a tunnel barrier (non-ideal coupling) described 
by a fixed scattering matrix ^i, say, the resulting S, constructed using the combination law 
of 5*0 and 5*1, is distributed according to PK. Brouwer's proof, being essentially a change of 
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variables from 5*0 to the final S, does not require stationarity, or ergodicity, or the maximum- 
entropy ansatz; however, it neglects evanescent modes between the barrier and the cavity. 
In other words, the reproducing property, which is fulfilled identically for the invariant 
measureiSiii, is not destroyed by the presence of the tunnel barriers. The latter certainly 
give rise to a nonzero {S), so that the direct processes described by this (S), being produced 
by the tunnel barriers, take place outside the cavity (see Fig. 2 in Ref. Il6| ). In contrast, 
when direct processes take place inside the system, it is not possible, in general, to write 
the total S as the combination of an 5*0 and a fixed Si, as required by Brouwer's analysis. 
Take, for instance, the system shown in Fig. Efd). If we had, say, a long "neck" between 
the small cavity and the big one, then we could define scattering matrices 5*1 for the former 
and 5*0 for the latter and combine them, disregarding evanescent modes, to obtain the total 
scattering matrix S. However, this is not the case for the actual system under study. As 
an approximation, we might think of assigning to the small and big cavities of the system 
of Fig. ini^d) the scattering matrices 5*1 and So, respectively, that would occur if we added 
the neck between the two; the total S obtained by combining these open-channel 5*1 and 5*0 
would represent an approximation to the actual problem; however, we are not in a position 
to know how close this approximation would be to the exact solution: we leave this open 
question for future investigation. Once again we seem to find that the valididty of PK for 
the systems studied in the previous section goes beyond the domain in which Brouwer's 
result was derived. 

Brouwer has also shown^^ that PK for the S matrix can be obtained from a Lorentzian 
ensemble of Hamiltonians with an arbitrary number of levels M. In the limit M oo the 
Lorentzian ensemble becomes equivalent to a Gaussian ensemble. In this limit, in which we 
believe that the Gaussian ensemble describes a chaotic cavity, the problem becomes once 
again stationary in energy. 

It thus seems that a derivation of PK -or at least of the reproducing property- for chaotic 
cavities with a general type of direct processes and in the absence of stationarity is, to our 
knowledge, still missing. 

When this work was completed, the present authors became aware of a study of the shot 
noise problem by D. Savin et al.,— and P. Braun et al.— in which results similar to those of 
our Sec. Ill G I have been obtained. 
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APPENDIX A: THE POLAR REPRESENTATION, THE INVARIANT MEA- 
SURE AND SOME STATISTICAL DISTRIBUTIONS FOR ONE AND TWO 
CHANNELS 

For completeness, we present the explicit parametrization of the S matrix in the polar 
representation for = 1 and N = 2 and some of its applications. 



1. The one-channel case, = 1. 





r 






t 





We write the two-dimensional S matrix in the polar representation as 
and the optical {S) as 

' {r) (f) 
it) {r') 

where the various entries are complex numbers. For (3=1 one has the restrictions 7 = 0; and 
5 = (3. The distribution of the conductance T = r for /5 = 1 can be reduced to quadratures. 



{S) 



(Al) 



(A2) 



with the result given in the text, Eq. p.llj) . 



a. The two-channel case, N = 2. 

The expressions given below are used in the present work when carrying out the numerical 
computations; since these are performed for the orthogonal case, /5 = 1, we restrict ourselves 
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to this universality class. For /3 = 1 we write the four- dimensional S matrix in the polar 
representation as 



S 





r 


t' 




t 





(A3) 



The reflection and transmission matrices r, t, etc., are two dimensional. The matrix r is 
two dimensional and diagonal: Tab = Ta^ab, with < < 1. The matrices v^^^ and w*-^-* are 
two-dimensional unitary matrices which can be written as 



yi _ fit) gi(««+7'')) ^/:^gi(aW+<5W) 



(5) 



The optical (5) is written as 

'(r) it') 
it) W) 

where the various entries are two-dimensional matrices. 
The invariant measure for f is given by 



(2t)' 



with the range of variation of the parameters 



< a^'\ f3^'\ 1^^, S^'^ < 271, 
< f('^ < 1 



(A4) 



(A5) 



(A6) 



(A7a) 
(A7b) 



and is normalized as J dn{v^^^) = 1. 

The joint probability distribution of Ti, T2 is given in Eq. ()2.19ap of the text: it is a 10- 
dimensional integral, with dn{v^^^) given in Eq. ()A6|) . the range of variation of the parameters 
being specified by flA7|) . 
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